Parity violation constraints using 2006-2007 QUaD CMB polarization spectra 
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We constrain parity- violating interactions to the surface of last scattering using spectra from the 
QUaD experiment's second and third seasons of observations by searching for a possible systematic 
rotation of the polarization directions of CMB photons. We measure the rotation angle due to 
such a possible "cosmological birefringence" to be 0.55° ± 0.82° (random) ±0.5° (systematic) using 
QUaD's 100 and 150 GHz TB and EB spectra over the multipole range 200 < t < 2000, consistent 
with null, and constrain Lorentz violating interactions to < 2 x 10 -43 GeV (68% confidence limit). 
This is the best constraint to date on electrodynamic parity violation on cosmological scales. 

PACS numbers: 11.30.Er, 98.80.Es, 98.70.Vc, 95.85.Bh, 95.30.Gv 



BACKGROUND 

Cosmic Microwave Background (CMB) polarization 
measurements at multipoles of £ > 20 are unaffected 
by reionization and are an effective means to probe for 
cosmological scale electrodynamic parity violation to the 
surface of last scattering. Using the CMB is particu- 
larly attractive because of the long path length to the 
surface of last scattering, the well-understood physics of 
the primordial universe that generated the CMB pho- 
tons, and two cross-spectra, the temperature-curl (TB) 
and gradient-curl (EB) cross-correlationSjthat should be 
null in a parity-conserving universe [2, 0, 0, 0, HI • As the 
effect should be frequency independent, measurements 
of the CMB at multiple frequencies can distinguish it 
from other EB correlation inducing effects like Faraday 
rotation from magnetic helds in the intergalactic medium 

The known parity violation in the weak force is suf- 
ficient motivation for investigating electrodynamic par- 
ity violation, but it has been shown that parity-violating 
interactions are a potential solution to the problem of 
baryon number asymmetry because they can be a signa- 
ture of CPT (charge-parity-time) violation in an expand- 
ing universe 



The effect arises by adding a Cherns-Simons term 
to the normal electrodynamic Lagrangian, violating 
Lorentz, P and CPT symmetries 1(J 11 1: 



C=- 7 F ltu F» v +p ll A v F>" 



(1) 



Here F^ u denotes the field tensor, F^ v is its dual, p M is 
an external vector, and A v the 4-vector potential. Non- 
zero time or space components of p^ induce a rotation of 
the polarization direction of each photon as it propagates 
from the surface of last scattering. This is equivalent to a 
local rotation of the Stokes parameters, Q and U, in the 
polarization maps made by CMB experiments, inducing 
gradient (E) to curl (B) mode mixing and therefore EB 
correlation. Lorentz violation can also be tested with 
these models [lfj, ll2| • In addition, models of quintessence 
can be probed by examining the EB and TB spectra for 



non-zero power 13 1 



QUaD was a 100 and 150 GHz bolometric polarimeter 
that made deep observations of the CMB from the South 
Pole during the austral winters of 2005 through 2007. A 
recent analysis of the second and third seasons of data 
from QUaD shows a series of acoustic peaks in the EE 
auto-spectra over the multipole range 200 < I < 2000 
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consistent with the A-CDM model of the universe 14 1. 
This dataset offers the strongest constraining power to 
date on cosmological scale parity-violating interactions. 
The QUaD collaboration maintains two code indepen- 
dent, but nearly algorithmically identical data analysis 
pipelines for the purposes of consistency checking. The 
results presented here use the 100 and 150 GHz spectra 
from the "alternative pipeline" described in section 6.8 
of Pryke et al. 14] for reasons of computational conve- 
nience, derived using a modified version of the MASTER 
CMB analysis method [HI . 



across our 23 reported bandpowers. 
£(£ + l)C[" 
fine the quantities below for each bandpower 



First, we assume 



l)C^ X ' obs is constant within a bandpower and de- 
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Debi 



Cj B ' obs cos (2Aa) - Cj E ' obs sin (2Aa) (8) 
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1 , n BB,obs 



Cf E ' ohs ) sin (4Aa) 



ANALYSIS 



Assuming that the CMB is a Gaussian random field, 
the entirety of its statistical properties can be described 
by the auto- and cross-correlation power spectra: 



We can then minimize x 2 (Aa) for the TB and EB 
combinations separately to estimate Aa : 



Cf Y — 2i+\ a fm a Jm (2) 

m 

where the ai m are the coefficients of the spherical har- 
monic decomposition of the temperature or polarization 
maps. X and Y here denote T, E or B for the respec- 
tive maps of temperature, gradient-polarization and curl- 
polarization modes. 

Normally the Cj B and Cf B are expected to be null 
because the spherical harmonic eigenfunctions Y^ and 
Y^ have parity (— 1)^ and has parity (— As- 
suming that there is a parity-violating effect in the elec- 
trodynamics equations that prefers one polarization to 
another over cosmological scales, let us denote the av- 
erage preferred rotation of the polarization direction of 
a photon from the surface of last scattering as it heads 
towards us as Aa. This corresponds to a rotation of the 
polarization directions in the maps [j], Q inducing E to 
B mixing, and therefore EB cross correlation. Likewise, 
since there is already TE cross correlation, TB cross cor- 
relation is also induced. Following Komatsu et al. (l5| . we 
assume that cosmological BB modes are zero to simplify 
the equations and maximize the likelihood of a detection: 



X 2 (Aa) = Y,Dtb,iM£D tb s (10) 

w 

X 2 (Aa) = J2 D EB,eM-}D EBi , (11) 

u 



We empirically measure the covariance matrix, M u > , of 
the bandpowers in each spectrum Debi and Dtb,£ from 
a set of simulated bandpowers combining realizations of 
A-CDM cosmology temperature and polarization fields 
for the signal component and accurate realizations of 
QUaD's instrumental noise. Our method utilizes a set of 
496 signal and noise Monte Carlo simulations from the 
analysis pipeline of QUaD. Pryke et al. [3| demonstrates 
the robustness of QUaD's simulation method against a 
variety of systematics tests. 
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sin(2Aa) 
Cf E cos 2 (2Aa) 
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(3) 
(4) 
(5) 
(6) 

(7) 



For the purposes of plotting and analysis, we can de- 
rive a theory-independent % 2 statistic to combine the 
first two and the last three equations separately to ob- 
tain an estimate of Aa, utilizing constraining power from 



FIG. 1: Aa measured from QUaD 150 GHz TB and EB spec- 
tra; histogram of simulations and red line for data. Histogram 
does not account for systematic error. Dotted line indicates 
total uncertainty assuming a Gaussian 0.5° systematic error. 



It is also possible to estimate Aa by measuring the quanti- 

or ,EB,abs n TB,obs 

on a per 



ties 



and 



bandpower basis, combining them using the covariances as mea- 
sured from simulations, and then applying inverse trigonometric 
functions. However, this is biased in the presence of noise. We 
thank an anonymous referee for suggesting our current method. 
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TABLE I: Column 1: Aa measurements from QUaD, including random and systematic errors. Column 2: Bias and standard 
errors on mean sampled from 496 signal-only simulations. Column 3: Column 2 subtracted from Column 1. Column 4: Scatter 
of signal-only simulations, indicating sample variance. Column 5: Fraction of signal+noise simulations where ||Aa|| exceeds 
that of data. 





Aa 


systematic 


bias-corrected Aa 


signal-only 


% simulations 


Spectrum 


(random and sys. errors) 


bias 


(random and sys. errors) 


simulation scatter 


exceeding 


150 GHz EB 


0.76°±0.92° ± 0.5° 


0.003°±0.003° 


0.76°±0.92° ±0.5° 


0.08° 


41.3% 


150 GHz TB 


1.19°±3.26° ± 0.5° 


0.025°±0.017° 


1.16°±3.26° ±0.5° 


0.37° 


71.5% 


100 GHz EB 


-3.74°±2.22° ± 0.5° 


0.011°±0.004° 


-3.75°±2.22° ±0.5° 


0.10° 


8.87% 


100 GHz TB 


3.72°±5.69° ± 0.5° 


0.073°±0.022° 


3.65°±5.69° ±0.5° 


0.50° 


52.2% 


150 GHz combined 


0.85°±0.94° ± 0.5° 


0.015°±0.003° 


0.83°±0.94° ±0.5° 


0.07° 


35.8% 


100 GHz combined 


-1.86°±2.24° ± 0.5° 


0.031°±0.005° 


-1.89°±2.24° ±0.5° 


0.11° 


38.7% 


100/150 combined 


0.56°±0.82° ± 0.5° 


0.011°±0.004° 


0.55°±0.82° ±0.5° 


0.08° 


49.6% 



150GHz EB parity violation limits 




150GHz TB parity violation limits 
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FIG. 2: 150 GHz Aa per bandpower derived from the EB 
spectrum. Note that in practice these points are combined 
before the final transformation to Aa — the purpose of this 
plot is to give a visual representation of the relative uncer- 
tainties across the bandpowers. 
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Figure Q] shows the results of this combination for the 
data (red line) and simulations (histogram) for 150GHz 



in both EB and TB. Overplotted is the total uncertainty 
assuming that the simulations reflect a normal distribu- 
tion and that the systematic error is 0.5°. It is clear 
that the observed data can easily be drawn from the set 
of simulations in which no parity-violating interactions 
have been included; we therefore conclude that there is 
no detection. 

To obtain a visual representation of a "Aa spectrum," 
We can also estimate the best fit for Aa on a per- 
bandpower basis by minimizing: 



(Aa) = J2 D TB,eM-}D TB/ (12) 

t 

X , 2 (Aa) - J2 D EB,iM~}D EB/ (13) 

The Aa spectrum using the EB, BB and EE spectra for 
150GHz is shown in figure [2] 



CURRENT LIMITS AND QUAD RESULTS 

Komatsu et. al. [HI report their limits from the 
WMAP 5 year high-^ data as Aa = -1.2° ± 2.2°. 
Other authors have found weak evidence for parity vi- 
olation by combining the WMAP 5 year data and data 
from the BOOMERanG balloon experiment, reporting 
Aa 



2.6° ±1.9° [llj. Carroll et. al. Carroll et al. 



derived constraints on Aa 10 high-redshift radio galaxies 
in 1990, yielding Aa = -0.6° ±1.5°. The best single 
redshift number, for 3C9 at z = 2.012, is Aa = 2° ± 3°. 

QUaD's results broken down by individual spectrum 
and frequency, as well as combined within and between 
frequencies, are shown in Table HI Reported errors are 
68.2% confidence limits as determined by the distribu- 
tion of signal and noise simulations. 150 GHz EB alone 
is significantly more constraining than any current result. 
At no frequency, nor in any spectrum, is there a signifi- 
cant detection. We also present values for Aa where the 
systematic bias induced by a combination of timestream 
filtering and the slightly different, non-aligned, and el- 
liptical nature of the beams of two orthogonally aligned 
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polarization sensitive detectors within a single feedhorn 
leading to temperature to polarization leakage has been 
quantified by signal-only simulations. This effect is dis- 
cussed in further detail in Hinderks et al. [13] • Note that 
in all frequencies and spectra this bias is an order of mag- 
nitude smaller than our random and systematic errors. 
After combination the EB spectra dominate the analysis 
and there is virtually no bias. These results are con- 
sistent with a constraint on isotropic Lorentz-violating 
interactions of fc ( ( ^ )00 < 2 x 1CT 43 GeV [H. 

SYSTEMATIC EFFECTS AND CHECKS 



that their contribution to the final uncertainty is small. 

As detailed in Table QJ our analysis of signal-only sim- 
ulations reveal a small bias in the recovered Aq values. 
In order to isolate the source of this bias, we have per- 
formed additional sets of signal-only simulations, includ- 
ing in isolation the effects of filtering, mis-aligned beams, 
uncertainties in detector alignment and cross-polar leak- 
age. The results from these tests confirm that a combi- 
nation of timestream filtering and beam mis-alignment is 
the source of the bias. Note that, although small com- 
pared to our noise-driven errors, our results do include a 
correction for the bias. 



The primary systematics concern is that there might 
be a systematic rotation of the true detector sensitivity 
angles, producing a false signal totally degenerate with 
that of parity violation; for example, a —3° systematic 
misalignment and a Act = — 3° true parity violation sig- 
nal would produce identical results. We have measured 
the overall rotation angle of our instrument using two 
methods. The first measures the polarization sensitiv- 
ity angle of each bolometer using a near field polariza- 
tion source. The second constrains the absolute angle 
of the focal plane by examining the measured offsets of 
the beams of each detector from the telescope pointing 
direction on an astronomical source. These two methods 
agree nearly exactly indicating that any systematic rota- 
tion of the bolometers within the focal plane structure is 
negligible. Given that there is no physical or mechanical 
reason to suspect such a rotation, and the uncertainties of 
the measurements, we conservatively assign a systematic 
uncertainty on the absolute rotation angle of the instru- 
ment of 0.5° and quote this value in the abstract and in 
Table H 

We have reanalyzed the entire dataset after inserting 
an artificial 2 degree local polarization rotation only in 
the data maps, resulting in a 2 degree shift after deriv- 
ing Aa identically to the procedure above, validating the 
analysis pipeline. 

A secondary concern is random scatter in the assumed 
detector angles. This is a different effect than a system- 
atic rotation of all of the detectors. The Monte Carlo 
simulation pipeline includes the injection of a degree of 
uncertainty about the true orientation of each polar- 
ization sensitive bolometer (PSB) into every simulation 
commensurate with the uncertainty of the measurements 
described in Hinderks et al. [17j |. Thus, when construct- 
ing a "fake focal plane" for signal-only simulations of a 
given CMB realization, we assign every bolometer a ran- 
dom deviation from its presumed angle at reconstruction, 
drawn from a Gaussian distribution with a = 1° . We also 
assume that the polarization grids are sensitive to the or- 
thogonal polarization direction at a level of 7 ± 3% and 
include this effect in the simulation pipeline. Signal-only 
simulations with and without these effects included show 



CONCLUSIONS 

We have presented the strongest constraints on par- 
ity violation to date. Assuming that there are no 
cosmological-scale parity violating interactions, we have 
also demonstrated that it is possible to understand the 
cumulative effects of detector misalignment uncertainties 
in polarization sensitive bolometer-based instruments to 
under 1° through a combination of analysis of primary 
CMB polarization data and lab measurements. This is 
of potential interest with respect to analysis of data from 
the High Frequency Instrument of the upcoming Planck 
Satellite. 
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